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INTRODUCTION 


The  Nuclea,' i Bunst  Detection  System  (NBDS)  net-ici  on  pnocessoig 
opticai  put  id  fin.  om  above  ground  nuctean  mapon  detonations  to  pno- 
vide  infionmation  to  battlefiield  com mandens.  The  system  penfionmance 
needs  to  be  venifiied  unden  nealistic  conditions  ofi  atmosphenic  tnans- 
panency,  cloud  nefitections,  and  tennain  shadowing.  A necentty  devel- 
oped concept  fion  a thenmo chemical  fitask  genenaton  appeans  to  ofifien 
pnomise  fion  a fiieldable  sounce.  The  pnimany  panpose  ofi  this  nepont 
is  to  evaluate  the  potential  fion  a low  cost,  neliable  MBPS  test  system. 

1.1  Rationale 

In  order  to  equip  battlefield  commanders  with  information 
on  the  times  and  locations  of  nuclear  weapon  detonations  in  wartime 
situations,  a Nuclear  Burst  Detection  System  has  been  developed  that 
relieson  (among  other  signals)  the  optical  pulse  signature  of  an  air 
burst  weapon.  This  report  is  concerned  solely  with  the  optical 
portion  of  the  NBDS,  and  the  desire  to  provide  an  optical  source  that 
will  trigger  and  test  the  system  over  a wide  range  of  realistic  con- 
ditions. Such  considerations  as  the  optical  pulse  risetime,  width 
and  intensity  observed  by  the  optical  detectors  in  an  NBDS  are  primarily 
functions  of  weapon  yield,  height  of  burst,  range  and  direction.  They 
will  be  modified,  perhaps  seriously,  by  atmospheric  transparency,  cloud 
shadowing  and/or  secondary  reflections,  terrain  shadowing  and  reflec- 
tions, and  the  effects  of  weather  (haze,  rain,  fog,  etc.).  Lightning 
discharges  (particularly  the  rare  super-lightning  type)  may  provide 
false  triggers.  A system  test  with  simulated  nuclear  pulses  would 
resolve  areas  of  uncertainty  in  NBDS  response.  It  is  the  purpose  of 
this  work  (drawing  on  recent  experimental  data  on  optical  flash 
generators)  to  calculate  the  theoretical  conditions  required  for  a 
valid  NBDS  test.  Furthermore,  test  system  configurations  of  likely 
use  in  NBDS  validation  are  recommended. 

1.2  Optical  Flash  Sources 

A source  for  the  NBDS  tests  requires  very  high  brightness 

in  the  pass  band  (visible  and  near  IR),  reasonably  wide  angular  diver- 
gence, and  controllable  risetime  and  duration.  These  conditions 
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may  rule  out  a number  of  conventional  choices  (e.g.  lasers).  Recently 
Science  Applications,  Inc.  (SAI)  has  completed  an  initial  development 
program  for  a large  scale  optical  flash  generator  based  on  high  temper- 
ature thermochemical  reactions.1,2  Tests  of  a field  proven  low  cost  unit 

g 

of  109  watts  peak  power  and  above  10  joules  of  optical  energy  have 
raised  the  possibility  of  use  of  this  device  for  NBDS  validation.  In 
this  report  an  assessment  of  the  potential  of  themochemi cal  reactions 
as  a source  of  optical  test  pulses  for  the  NBDS  is  made. 

Development  of  the  thermochemical  optical  flash  generator  prior 

to  this  work  concentrated  on  the  aspects  of  one  specific  reaction,  that 

of  aluminum  oxidation  at  high  temperatures.  Initial  theoretical  research 

had  isolated  metallic  oxidization  as  a potential  source  of  large  scale 

optical  flash  output.  The  temperature  of  such  a radiating  system  is 

limited  on  the  upper  side  by  the  vaporization  temperature  of  the  metallic 

oxide.  Alumina  (Al^O^)  with  a vaporization  temperature  of  around  3800°K 

(see  Appendix)  would  allow  radiation  fluxes  of  about  200  cal/cm2sec  if  a 

sufficiently  optically  thick  cloud  of  A1  and  0^  can  be  created  and  ignited 

at  high  temperature  forming  A1 2^3  near  the  vaporization  point.  Aluminum 

is  readily  available  in  various  finely  divided  forms, is  inexpensive  and  is 

non-toxic.  Thus,  previous  work  proceeded  to  ascertain  the  experimental 

details  of  a practical  generator.  A series  of  about  50  tests  were  performed 

in  1976  and  reported  by  Cockayne  et  al.1  Later  a series  of  over  100 

thermochemical  reaction  test  firings  were  conducted  in  a continuing  develop- 
9 

ment  series  . These  and  open  literature  sources  on  thermochemical 
reaction  tests  form  a basis  for  theoretical  development  of  the  application 
to  NBDS  validation. 

1.3  Approach 

The  usefulness  of  any  optical  source  for  NBDS  testing  depends 
on  the  trigger requi rements  set  by  the  sensor  system  and  on  the  objec- 
tives set  for  system  validation.  The  approach  in  this  report  will  be 
to  review  and  relate  the  NBDS  requirements  for  the  optical  output  of 
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both  nuclear  and  thermochemical  devices.  To  do  this,  development  of  the 
expected  output  of  large  scale  thermochemical  reactions  will  be  inferred 
as  scaling  laws  based  on  data. 


In  summary  the  key  issues  to  resolve  are  these: 

1)  Does  the  thermochemical  generator  provide  a 
trigger  for  the  NBDS? 

Z)  What  are  the  basic  equations  that  describe 
the  thermochemical  generator  output? 

3)  Are  there  alternate  schemes  of  calibration 
that  holu  promise? 

The  purpose  of  this  report  is  to  answer  these  questions. 


fMrw. 


II 


OPTICAL  OUTPUT  OF  NUCLEAR  WEAPONS  AND  NBDS  RESPONSE  PARAMETERS 


The  optica?  output  ofi  cunbutit  nucteaA  weapon*  watt  be  pcuuz- 
metaized  veuui  yield,  xange  and  height  o&  bun*t,  The*e  pnedlation* 
will  be  contra* ted  with  the  expected  Viiggeh  aeqwiacment*  o{,  the  MBPS 
ii/j  tem. 

II. 1 Nuclear  Weapon  Optical  Output 

The  first  task  of  the  reported  work  requires  a characteriza- 
tion of  the  optical  (thermal)  waveforms  of  nuclear  bursts.  It  is 
well  known  that  airburst  nuclear  weapons  have  a unique  optical 
signature  consisting  of  a first  rapid  pulse  followed  by  a slower 
and  much  longer  pulse.  The  time  of  the  peak  of  the  first  pulse  is 

t,  > similarly  t0m,  is  the  time  of  the  second  peak,  and  t . is  the 
I max  2max  min 

time  from  initiation  to  the  minimum  between  pulses.  The  most  general 

review  of  nuclear  phenomonology  - that  of  Glasstone  - quotes  the  follow- 

3 

ing  relationships  for  the  thermal  output: 


II  .A. 

P 

max 

4W 2 cal/sec 

II.  B. 

^tot 

W/3 

II  .C. 

t . = 

min 

2 . 5W 2 msec 

II  .D. 

t2max 

= 32W 2 msec 

12 

where  W is  the  weapon  yield  in  KT  (1KT  = 10  cal).  The  first  pulse 
contains  approximately  1 percent  of  the  total  thermal  output.  Figure  1 
is  taken  from  Reference  3 and  contains  details  of  the  normalized  pulse 
shape  (second  maximum  only)  and  total  output  for  various  yields. 


In  order  to  calculate  flux  (Q)  and  fluence  (Q)  at  a distant 
point,  ground  range  R from  a burst  of  H height,  we  need  additional 

formulae*  fhpv  a rp • 


Q = 

W 

1 o8 

(Cal/cm2) 

1 2tt 

(H2+R2) 

Q = 

U 

108 

(Cal/cm2  sec) 

TT 

( _i_n2  \ 

No  atmospheric  transmission  loss  is  included  in  these  formulae. 

Additional  estimates  of  the  time  spectrum  have  been  obtained  from 
Reference  4.  Figure  2 shows  an  example  of  one  of  these  and.  Table  1, 
a compilation  of  pertinent  data  at  several  yields.  Figures  3 and  4 
summarize  the  tmin  and  t2mdx  calculation  (Table  1)  and  formula  of  Ref- 
erence 3. 

1 1. 2 NBDS  Operation 

Optical  pulses  are  detected  by  a Si  based  detector  whose  response 
function  is  shown  in  the  accompanying  figure  (Figure  5).  There  are  four 
real  time  hardware  discriminants  that  determine  whether  an  event  is 
likely  to  be  a nuclear  event  and  which,  therefore,  trigger  the  retention 
of  0.32  msec  of  data  in  the  NBDS  memory , accumulation  of  additional  data, 
and  activation  of  the  microprocessor  discriminant  analysis.  These 
discriminants  are  listed  below. 

Real  Time  Discriminants 

1)  VHF  Negative  Pulse  exceeding  a specific  threshold  5 v/m 

2)  VHF  Arrival  within  1 millisecond  of  optical  trigger 

3)  Optical  signal  trigger  level  >6.3  x 10  watts/cm 

4)  Optical  rate  of  rise  >2.1  x 10~3  watts/cm^/sec 

The  first  two  of  the  real  time  discriminants  are  electromagnetic 
pulse  discriminants  which  are  not  relevant  to  our  problem.  The  last 
two  merely  indicate  a significant  increase  above  an  ambient  light  thres- 
hold and  a rate  of  rise  consistent  with,  but  not  exclusively  associated 
with,  a nuclear  event. 
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Figure  2.  In-Band  Power  at  Source  for  Two  Broad  Band  Silicon  Sensors 


Relative  Response 


Microprocessor  Discriminants 


1) 

2) 

3) 

4) 
*) 
6) 

7) 

8) 
9) 

10) 

11) 

12) 


First  sample 


tl.nax  3-6  ",sec 


6.3  x 10  watts/cnr 


time  of  (first  max/6. 25)>  0.5  msec 


5 t 


lmax 


min 


70  t 


lmax 


4 t 
t 


t. 


'min  2max 

“min  199  insec 
Width  of  min  < 8 t 


20  t 


mi  n 


min 


W6-25 


88  of  the  264  samples  before  100  msec 
Flat  long  duration  first  peak  pulses  are  rejected 
First  max  > avg.  amplitude  of  up  to  128  samples 
Second  max  > 6.25  ■ min  amplitude 

First  ten  samples  after  second  max  > second  max  : (6.25/3) 


The  first  of  the  mi croprocessor  discriminants  is  presently 
redundant  but  could  be  varied  to  exclude  low  yield,  distant  bursts 
from  further  consideration.  The  second  is  a check  on  the  rise  time 
of  the  first  maximum  which  is  applicable  to  all  normal  yields.  The 
third  is  actually  performed  by  checking  that  the  amplitude  of  the  pulse 
at  0.5  msec  is  greater  than  1/6.25  of  the  first  maximum.  (The  factor 
6.25  is  equivalent  to  the  difference  between  any  consecutive  32  analog 
to  digital  conversion  steps  in  the  sensitivity  range  of  the  instrument 
and  is  used  by  the  microprocessor  as  a consistent  figure  for  measuring 
amplitude  variations.)  If  this  is  the  case  then  the  first  maximum  has 
not  formed  by  0.5  msec  and  consequently  the  first  maximum  is  broader 
than  this  value.  The  fourth  and  fifth  determine  the  arrival  times  of 
the  minimum  and  second  maximum  in  terms,  ultimately,  of  the  width  of 
the  first  maximum  and  are  consistent  with  all  yields  within  the  range 
of  0.1  KT  to  1 MT.  The  sixth  is  satisfied  by  all  normal  yields  and  is 
an  absolute  check  on  the  arrival  time  of  the  minimum  rather  than  one 
related  to  the  first  peak  arrival  time  ( t jmax ) as  in  the  fourth.  The 
seventh  and  eight  are  related  to  the  shape  of  the  first  pulse.  The 
eighth  says  specifically  that  a third  of  the  registers  covering  the 
first  100  msec  (not  a third  of  the  time  because  the  register  time 
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interval  is  not  constant)  have  amplitudes  within  the  first  thirty-two 
range  sensitivity  steps  below  the  first  maximum  value  and  serves  as  a 
check  that  the  first  peak  is  not  too  irregular  in  shape.  The  seventh 
specifically  checks  that  fewer  than  36  of  the  registers  before  100  msec 
are  within  any  one  sensitivity  range  which  would  make  the  pulse  too  flat. 
The  ninth  specifies  that  the  second  peak  must  start  before  1800  msec 
(time  sample  #396)  whereas  discriminants  5 and  6 would  specify  that 
the  second  peak  must  occur  before  2000  msec.  The  ninth  would  therefore 
appear  to  be  redundant  at  present  but,  as  with  the  first,  could  be 
modified  to  increase  the  NBDS  discriminant  capability  in  the  future. 

While  the  eighth  insures  that  the  first  pulse  had  a significant  number 
of  registers  with  an  acceptable  amplitude,  the  tenth  checks  that  the 
pulse  is  also  not  too  uniform  or  flat.  The  eleventh  and  twelth  are 
designed  to  exclude  lightning  flashes  which  have  a smaller  second 
maximum  in  relation  to  the  first  maximum  than  a nuclear  burst  and, 
due  to  the  elongated  geometry  of  the  source,  cool  rapidly  to  provide 
a narrow  second  pulse. 

The  determination  of  whether  or  not  the  burst  is  a ground  burst 
is  based  on  the  definition  of  a ground  burst  as  one  whose  fireball 
intercepts  the  surface  of  the  earth.  Such  an  occurrence  has  at  least 
two  effects.  First,  the  fireball  will  cool  more  rapidly  as  it  heats 
the  ground  and  second,  the  fireball  will  never  reach  a second  maximum 
as  high  as  it  otherwise  would.  These  two  points  are  addressed  by  the 
Ground  Burst  Discriminants 


» l2n,ax  i 9-5916  x ‘min 
and/or 

2)  Second  max  < 1.8018  x First  max 
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III  OPTICAL  PULSE  SIMULATION  TECHNIQUES 

Simulation  technique b bait’d  on  thermochemical  generatorb 
uyi.lt  be  dlbcubbed.  Theoretical  argument b fior  the  gene-rat  beating 
laid  faound  to  hold  in  the  reactlonb  studied  to  date  bhow  that  the 
output  and  tame  duration  oh  the  pulbe  can  be  varied  within  llmitb. 
Contrabtb  with  other  bourceb  show  that  hlabh  tempi  albo  hold 
protnibe  hor  providing  at  leobt  part  oh  the  required  optical  big- 
nature.  Optical  enhancement  oh  the  output  oh  any  omnidirectional 
radiating  bource  may  prove  to  be  cat.  elective  in  hie-td  tebtb. 
Detailed  uuderb  tanding  oh  the  internal  mechanaimi  oh  h-ireball  develop- 
ment promlbei  additional  wayb  to  tailor  the  output  pulbe  oh  the  gen- 
erator. 


I II. I General  Background 

No  device  exists  at  present  that  has  been  used  success- 
fully to  simulate  nuclear  pulses  to  NBDS  detectors,  although 
accidental  triggering  of  NBDS  due  to  naturally  occurring  "super- 
lightning" has  been  reported.  The  use  of  a thermochemical  generator 
based  on  a metallic  oxidization  reaction  has  been  proposed  as  a 
calibration  technique.  A brief  review  of  its  operation  will  be 
presented.  Consider  the  metallic  oxidization  reaction: 


M + aO. 


M0„  + heat 

2a 


The  highest  final  temperature  practically  obtainable  due 
to  this  reaction  is  that  of  the  vaporization  temperature  (T  ) of 
the  oxide  MO^,  since  above  this  temperature  latent  heat  effects 
would  dominate.  Experience  has  shown  that  if  the  reaction  can  be 
initiated  either  at  Tf  or  the  lower  melting  point  temperature  it  will 
self-sustain  at  T . A thermochemical  flash  generator  is  most  radiation 
efficient  when  there  is  a large  surface  to  volume  ratio,  a high  tempera- 
ture T , and  a sufficient  optical  density  of  the  reaction  products.  Thus 
a system  of  metallic  dust  or  fine  wires  mixed  with  oxygen,  and  ignited 
simultaneously  at  many  locations  is  preferred. 

The  initial  energy  is  supplied  in  chemical  form  by  the 
fuel  (M  and  0^).  The  optical  fluence  is  then  proportional  to  W, 
the  weight  of  M, if  sufficient  0^  for  reaction  is  assumed.  The 
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optical  flux  per  unit  area  of  the  reacting  system  (a  fireball)  is 
governed  by  the  blackbody  equation  if  the  system  is  optically  thick. 

4* 

Thus  total  radiant  heat  loss  is  proportional  to  T and  area.  If  T 
is  constrained  to  T , then  the  total  flux  of  radiation,  Q,  is  proportional 
to  the  square  of  the  system  radius.  If  the  fuel  consumption  rate  exceeds 
the  ability  of  the  system  to  radiate,  the  system  will  expands losing 
energy  proportional  to  the  time  rate  of  change  of  the  radius.  When 
radiation  equilibrium  is  established  the  volumetric  expansion  will 
cease  although  the  apparent  size  may  increase  due  to  buoyancy  effects 
and  local  turbulence  at  the  fireball  edge. 

Since  the  energy  loss  in  expansion  is  proportional  to  volume 

it  can  be  stated  that  the  stabilized  fireball  radius  will  be  pro- 
1/3 

portional  to  W . Consequently,  since  Q is  proportional  to  W,  and 
. 2/3 

Q to  area  W ),  the  characteristic  time  <t>  of  the  optical  output 

1 /3 

is  proportional  to  W , assuming  a constant  radiation  temperature, 
optically  dense  reaction  products,  and  ignition  and  burning  times 
less  than  the  characteristic  time. 


Two  distinct  types  of  thermochemical  generation  schemes  are 
possible.  The  entire  amount  of  fuel  (M  + 02)  may  be  ignited  at  once, 
or  M may  be  metered  into  the  02  volume  after  initial  ignition.  Recall 
that  M is  in  the  form  of  a fine  powder,  facilitating  its  injection. 

In  the  latter  case  the  additional  fuel  will  have  varying  effects, 
instantaneous  flux  will  be  lower  due  to  the  absorption  of  energy 
in  raising  the  fuel  to  ignition  temperature;  however,  the  reaction 
time  may  be  lengthened  arbitrarily,  and  the  expansion  of  the  fireball 
controlled  by  the  proper  metering  of  fuel. 

In  summary  the  key  considerations  for  a thermochemical 
generator  are: 

• Choice  of  the  metallic  system  M02a  on  the  basis  of  evolved 
heat  and  vaporization  temperature 

• Choice  of  the  scale  of  the  reaction  W 

• Choice  of  the  metering  rate  of  fuel,  including  ignition 
of  all  fuel  at  initiation. 


*A  silicon  dector  response  cutoff  at  1.1  pm  causes  a loss  of  the  infra- 
red tail  so  that  T5  is  a better  estimate  of  the  detected  or  passed  power. 
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The  investigation  of  the  reaction  system  will  follow  these  considerations. 

A review  of  selected  thermochemical  reactants  was  made  in  Reference 
1 and  is  reported  in  Table  2. 


III. 2 Comparison  of  Thermochemical  Pulse  With  Nuclear  Pulse 

Reference  1 reports  on  the  current  status  of  the  SAI  developed 
thermochemical  generator.  The  following  are  typical  parameters  for  the 
largest  unit  fielded  as  of  March  1977: 


Pmav  = 10^  watts  = 2.4  x 103  cal/sec 

Ilia  a 


Etotal  = 1'1  x 10  cal 


tm,v  = 10  to  50  msec 

Ilia  A 


III. A 


III  .B 


III  .C 


Assuming  an  isotropic  radiation  distribution  the  sensor  must  be  at  a 
range  of  no  greater  than  224  km  to  satisfy  microprocessor  discriminant 
number  1.  The  "generator"  has  no  first  maximum,  of  course;  and,  thus, 
the  other  criteria  cannot  be  met  at  present  with  a single  unit.  With  a 
single  reaction  the  pulse  shape  of  a typical  unit  will  simulate  that  of 
the  second  maximum  of  a small  yield  device,  or  that  of  the  first  maximum 
of  a large  yield  device. 

Leaving  aside  questions  of  the  double  pulse  device,  we  wish  to 
estimate  the  "nuclear  equivalence"  of  the  pulse.  First,  let  us  assume 
the  pulse  is  the  main  or  second  pulse.  A typical  thermochemical  pulse 
is  of  10  to  50  msec  duration  to  peak,  with  the  shape  controllable  to 
some  extent.  Inverting  the  familiar  relationship 


Snax  = 32  W(KT)  msec 


III.D 
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Table  2 


Energy  Release  of  Selected  Thermochemical 
Reaction  Mixtures 


HEAT  OF 

REACTION 

APPROXIMATE 

FLAME 

TEMPERATURE 

COMPOSITION 

PARTICLE 

kcal/g 

cal /cc02 

(°K) 

Sol id-Gas 

A1  + 02 

A12°3 

7. 4a 

11.1 

3900 

Zr  + 02 

Zr02 

2.9a 

10.8 

4400 

Mg  + 02 

MgO 

5.9a 

12.0 

3100 

Fe  + 02 

FeO,  Fe^O^ 

1 . 6a 

5.6 

- 

Ti  + 02 

Ti  02 

4.7a 

9.4 

- 

Be  + 02 

BeO 

15. 93 

11.9 

- 

B + 02 

B2°3 

74. 0a 

8.4 

- 

Gaseous 

A1  (CH3)3  + 02 

CO 

o 

C\J 

< 

- 

4.6 

- 

Propane  + 02 

- 

- 

3.3 

- 

Sol  id 

A1  + LiCA04 

ai2o3 

3.0 

7 ,655b 

3900 

aBased  on  weight  of  metal  fuel. 
BSolid  mixture,  cal/cc. 
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and  solving  for  W yields: 


t (msec) 

W(KT)  = III.E 

W = 0.1  to  2.4  KT 

for  the  range  of  t „ specified  by  III.C. 

lila  X 

Another  method  of  calculating  the  yield  would  be  to  note  that, 
12 

of  the  10  cal/KT  released  by  a nuclear  device,  a third  of  the  energy 
is  typically  thermal  and  most  of  this  (.0.95)  appears  in  the  second 
pulse.  Then  III.B  leads  to: 

W = 3 x 10'5KT. 

Finally,  the  relation  between  the  peak  flux  for  the  second 
maximum  and  fluence  is: 

p«xl4<ll|l!  iii.f 

and  using  III. A we  find 

W £ | x 10"8  KT 

Thus,  the  pulse  width,  peak  power,  and  total  energy  of  the 

optical  output  of  the  simulator  are  not  consistent  yield  indicators 

for  a nuclear  weapon.  Yield  equivalents,  depending  on  method  of  calcu- 

-9 

lation, range  from  2 x 10  KT  to  2.4  KT.  Methods  of  enhancement  that 
direct  more  flux  toward  the  sensor  than  that  obtained  from  a radiating 
fireball  would  be  required  to  bring  these  values  into  synchronism. 
Additional  tailoring  of  the  pulse  shape  may  also  be  required  to  simul- 
taneously match  the  constraints  if  such  consistency  is  to  be  demanded 
by  an  NBDS  system. 

On  the  other  hand,  if  we  assume  the  simulator  is  employed  for 
the  first  maximum  then: 


III  .G 

min  * 2-5  * 10"3  “1/2 

III  .H 

lmax  = 2 * 10'4  “1/4 

III  .J 

and  are  appropriate  for  use  with  the  information  in  III. A and  III.C. 
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We  then  obtain 

W%  Jx  10'4  KT 
W ^ 2500  KT 

O- 

and  W % 104  KT 

% 

respectively.  And  we  once  again  find  a range  of  values  for  W which 

Q 

covers  almost  a factor  of  10  . 

The  overall  conclusion  is  that  a source  for  a supplemental 
pulse  is  needed.  This  may  be  a fast  rising  pulse  to  simulate  first 
maximum,  or  a very  long  high  intensity  source  to  simulate  second 
maximum.  The  supplemental  pulse  may  be  from  a revised  version  of  the 
thermochemical  reaction  system  or  from  a different  source. 

III. 3 Review  of  Potential  Optical  Sources  for  the  Second  Pulse 

An  initial  review  of  potential  sources  for  the  supplemental 
pulse  has  been  conducted.  The  following  candidates  were  considered: 

Black  Body 

a)  Thermochemical 

• ai/o2 

• Zr/02 

• Others 

b)  Shock  Radiation 

• Cylindrical  - in  02  or  Air 

- in  Argon 


23 


• Supercri tical  shocks  by  seeding  a gas  with  an 
optically  thick  medium 

• Multiple  interacting  shocks 

• Implosions 


Quasi  Block  Body 
Flash  Lamps  - Xenon 

• Pulsed 

• CW 

Line  Radiation 
Lasers 

• Pulsed 

• Modulated 
Hybrids 


For  the  thermochemical  generator  to  be  successfully  modified,  theoretical 
investigation  of  its  performance  is  required;  this  is  the  main  subject 
of  this  effort.  A variety  of  metals  may  be  considered.  Zirconium  has 
a higher  flame  temperature  and,  thus,  may  produce  a brighter  output 
(see  Reference  1,  page  9 et  seg.).  Shock  radiation,  possibly  a by- 
product of  the  ignition  sequence, may  be  a potential  source.  Seeding  of 
the  shocked  medium  to  increase  radiant  output  is  possible. 

Flash  lamps  seem  to  offer  the  next  great  promise  in  terms  of 
pulse  shaping  and  brightness.  They  were  selected  for  further  investi- 
gation as  a back-up  to  the  thermochemical  system.  Lasers  seem  unsuit- 
able in  comparison. 

The  sources  listed  above  differ  in  source  characterization, 
size  and  brightness.  To  assess  the  desirability  of  potential  sources, 
some  considerations  of  optical  transmission  to  a distant  detector  were 
made  for  two  cases,  a hot  radiating  ball  and  a small  flash  source 
with  optical  enhancement  in  the  propagation  direction. 
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In  the  first  case: 


p oc  j^(-) 

detector  \R  / 

when  T is  the  temperature  of  an  assumed  black  body  ball  of  radius  r, 
radiating  isotropically  to  a detector  at  R. 

The  second  case  is  approximated  by  a flash  source  of  size 
r , at  the  focus  of  an  »f  = 1 parabolic  mirror  of  radius  r'.  The 
object  size  is  r'  + R sinot,  when  a % r /r'.  For  large  R this  reduces 
to 

p cc  (—  ) 

detector  \R  / 

If  we  assume  that  the  size  of  the  first  fireball  is  equal  to  that  of 
an  #f  = 1 parabola  then  the  relative  advantage  is  the  ratio  of  the 
source  temperatures  to  the  fourth  power.  Since  structural  supports  are 
built  in  the  first  case  to  approximately  the  fireball  dimensions,  and  in 
the  second  to  the  size  of  the  parabola,  an  assumption  of  equal  costs  is 
justified.  A much  more  detailed  cost/benefit  analysis  will  have  to  be 
made  eventually,  but  this  first  look  indicates  that  a small  thermo- 
chemical source,  a flash  lamp  or  other  small  source  enhanced  with  optics 
may  be  potentially  useful.  A second  result  of  this  calculation  is  to 
indicate  the  elative  worth  of  optical  enhancement  with  any  source. 

1 1 1. 4 Details  of  Flash  Lamp  Systems 

A brief  review  of  flash  lamp  systems  was  made  in  order  to  obtain 
comparative  information  for  analysis  of  cost  benefit.  Typical  linear 
or  helical  quartz  tube  flash  lamp  characteristics  depend  mainly  on 
excitation  energy.  Designs  of  up  to  lm  in  length  and  2 cm  in  diameter 
are  standard.  Typical  noble  gas  filling  pressures  range  between  0.2 
and  1.0  atmosphere,  with  Xenon,  Krypton  and  Argon  being  the  most  commonly 
used  elements.  Water  or  gas  cooling  may  be  provided,  or  natural  air 
cooling  used  if  repetition  rates  are  not  severe.  Typical  parameters 
for  large  (lm  long)  flashlamps  are  color  temperature  up  to  9400°K,  high 
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emissivity,  particularly  at  longer  wavelengths,  and  costs  of  about 
$300  to  $600  in  quantity  purchase.  Such  lamps  are  approximately 
50%  optically  efficient,  based  on  electrical  input,  and  may  be 
cycled  many  thousands  of  times  if  not  run  near  peak  output. 

Rough  cost  estimates  indicate  about  $1  OK  per  module  (in  quantity 
purchase),  including  power  support  for  a 4000  cal/pulse  optical 
output. 


IV 


COMPARISON  OF  THEORY  TO  TEST  DATA  - THERMOCHEMICAL  GENERATOR 


The  general  theoretical  arguments  (5 on.  fireball  Ignition  and 
development  are  compared  with  three  sets  oft  test  data  on  thermo  chemical 
reaction  te sts.  The  (\irst  series  o(>  fairly  large,  explosively  ignited 
tests  show  good  agreement  with  the  theoretical  predictions . A second 
set  o ^ small  scale  test  data  is  then  used  to  confirm  this  over  a wider 
range  ofi  parameters. 


IV. 1 Test  Data  Taken  Prior  to  March,  1977 

In  this  section  we  review  the  salient  test  data  of  Reference  1 
for  comparison  to  theoretical  arguments  developed  above.  We  discuss 
initially  the  comparison  with  the  general  conclusions,  then  refer  to 
quantitative  comparison  of  theory  with  experiment. 

From  the  general  conclusions  (Section  5,  Ref.  1)  we  note: 

a)  Fluence  scales  linearly  with  the  weight  of  aluminum, 

(all  other  parameter  ratios  being  held  constant)  in 
agreement  with  the  basic  equation  for  the  device. 


b)  Peak  flux  scales  linearly  with  the  weight  of  aluminum 

implying  small  variation  in  the  fundamental  output 

pulse  width  over  the  range  of  aluminum  weights  studied. 

Since  the  larger  tests  were  made  up  of  smaller  models 

fired  simultaneously,  this  is  not  surprising.  Later 

work  over  a wider  range  of  aluminum  weights  more  clearly 
1 /3 

revealed  the  W dependency  of  pulse  width. 


c)  Shock  effect  proportional  to  output  was  noted.  The 
shock  source  is  postulated  to  be  the  primacord  igniter. 
Since  this  is  proportional  to  output  in  the  modular 
scheme  used,  the  results  follow. 

1 73 

d)  Fireball  radius  is  found  to  be  proportional  to  W , 
where  W is  the  weight  of  aluminum  powder.  This  is 
consistent  with  a fraction  of  the  fireball  energy 
being  dissipated  in  fireball  expansion.  This  fact 
is  corroborated  by  the  relative  radiative  efficiency 


r 
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(^10  ).  Modularity  of  the  test  devices  also  would 
tend  to  produce  the  effect. 


e)  The  peak  fireball  temperature  of  about  3600°K  is 

consistent  with  that  of  the  vaporization  temperature 
of  alumina  (Al^O^)  while  the  temperature  at  fireball 
expansion  stabilization  (2400°K)  is  consistent  with 
the  alumina  freezing  point. 

Analysis  of  fireball  data  in  Reference  1 evolved  three  empirical 
scaling  laws,  for  fireball  radius  at  stabilization, 

rf(m)  = 2.5  W(Kg)  1/3 

for  the  time  of  fireball  stabilization, 

tf (sec)  = W ( Kg ) 1/3/16 

and  for  fireball  radius  as  a function  of  time  to  stabilization 
r(m)  = 5 (W(Kg)  t(sec))3^ 

An  average  optical  efficiency  for  the  fireball  of  10^  cal/gmAl  was 
measured.  Eight  shots  of  sizes  from  0.5  gm  to  nearly  15  Kg  were  used  to 
calculate  these  values. 

IV. 2 Test  Data  on  Small  Scale  Shots 

Over  100  small  scale  shots  of  thermochemical  reaction  systems 

2 

using  a variety  of  configurations  have  been  reported.  Since  the  optical 
output  was  detected  with  a single  photodiode  and  not  corrected  for  fireball 
size  or  temperature  effects, only  a very  limited  amount  of  information  can 
be  drawn  from  these  tests.  Since  these  tests  and  the  previous  series 
cover  a fairly  wide  range  of  aluminum  weight,  investigation  of  optical 
pulse  duration  scaling  can  be  made.  In  addition, since  the  upper  limit 
of  particle  size  was  specified,  search  for  particle  size  effects  may  be 
made.  However,  since  the  particle  size  distribution  is  unknown,  and 
since  several  different  grades  of  powder  were  employed,  no  firm  conclu- 
sions can  be  drawn.  In  any  case,  about  40  small  scale  tests  using 
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explosive  ignition  were  studied  for  pulse  length  effects.  Of  these, 

about  1/4  (11/41)  were  duds  and  were  ignored  in  the  analysis.  The 

average  pulse  widths  for  the  remaining  tests  are  plotted  on  Figure  6 

where  data  from  6,  18  and  60  micron  upper  limit  (of  5 gm  aluminum 

weight)  particle  diameter  tests  are  indicated  separately  along  with  a 

typical  error.  In  addi tion, tests  at  100  gm  weight  are  indicated 

individually.  Pulse  widths  from  the  larger  shots  performed  in  the 

first  test  sequence  are  plotted  also.  Despite  considerable  data  scatter 
1/3 

the  general  W trend  seems  well  verified.  The  data  scatter  may  be 
due  to  the  variety  of  aluminum  powders  and  other  experimental  variables 
or  to  the  lack  of  instrument  corrections. 

Investigation  of  the  balance  between  the  fundamental  expansion 
time  of  the  fireball  with  the  ignition  and  burning  times  of  the  particles, 
and  of  the  optical  thickness  of  the  fireball  indicate  allowed  regions  of 
reaction  parameters.  As  an  example  of  this  analysis  consider  the  interplay 
of  particle  ignition  plus  burn  with  fireball  expansion  to  peak  output 
(at  3000°K).  From  Appendix  A we  derive  on  general  grounds  (Eq.  L) 

t (sec)  = 0.0063  (W(kg ) ) 1/3 

and  that  the  particle  ignition  time  is  given  by,  Eq.  U, 

t.(sec)  = 2.5  x 10"6  (a(u))2 
while  burning  time  is  given  as,  Eq.  V, 

tD  (sec)  = 0.0002  a(p)  - .0004 

D 

Equating  the  numerical  value  of  the  expressions  for  ignition  plus 
burning  with  expansion  to  peak  output  yields  the  constraint  equation 

a(p)  = (2.52  x 103  (W(kg))1/3  + 1760) 1/2  - 40 

This  formula  is  evaluated  in  Table  3. 


Table  3 


Evaluation  of  the  Constraint  that  Particles 
Must  Ignite  and  Burn  in  the  Characteristic 
Expansion  Time  of  the  Fireball 

i-Ug)  a(u) 


0 

2.0 

0.0001 

3.3 

0.001 

4.9 

0.01 

8.0 

0.1 

14.0 

1 .0 

25.0 

10.0 

45.0 

100.00 

76.0 

Likewise  another  constraint  equation  can  be  derived  by  requiring  that 
the  fireball  radius  at  stabilization  be  one  mean  free  path  thick  or: 


Using  equations  derived  in  the  Appendix  we  obtain 
a(p)  = 23(W(kg) ) 1/3 

These  two  constraints  are  indicated  on  Figure  7 along  with  specific 
points  indicating  parametric  location  of  tests  conducted  to  date. 

Tests  far  from  the  criteria  tended  to  have  lower  efficiency  than  those 
on  or  near  the  criteria.  Tests  where  the  characteristic  time  for 
fireball  expansion  is  much  less  than  that  required  for  ignition  and 
burning  tended  to  be  less  efficient  than  one  might  expect.  In  the 
extreme  opposite  case  the  tests  also  show  reduced  radiation  efficiency. 
This  may  be  due  to  the  fact  that  more  rapid  ignition  and  burning  drives 
the  expansion  and  decreases  radiation  efficiency. 

IV. 3 Additional  Test  Data 

The  most  recent  series  of  tests  performed  by  J.  Dishon  (Refer- 
ence 2)  have  consisted  mostly  of  shots  using  a configuration  of  pre-mixed 
aluminum  with  oxygen.  This  set-up  is  also  capable  of  injecting  addi- 
tional aluminum  into  an  ignited  fireball. 
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In  the  case  of  the  pre-mixed  system  an  optically  thin  mixture 
is  desirable  for  longer  pulse  duration  and  the  opposite  is  true  for  a 
shorter  pulse.  In  the  case  of  an  optically  thinner  mixture,  higher 
optical  efficiency  is  anticipated  due  to  the  more  effective  radiation 
from  within.  In  optically  thicker  mixtures  expansion  competes  with 
radiation,  lowering  optical  efficiency.  Table  4 summarizes  the  quali- 
tative features  of  several  test  configurations  and  indicates  the 
expected  scaling  laws.  References  5 and  6 summarize  NBDS  data  obtained 
in  these  tests. 

Table  5 reports  scaled  fluence  as  a function  of  oxygen  to 
aluminum  ratio  and  particle  size  for  metered  flow  conditions  (System  C 
of  reference  2).  This  illustrates  the  data  trends  to  decreasing 
efficiency  with  particle  size.  These  data  indicate  averaqe  fluences  of 
2000  to  3000  cal/gm  are  attainable. 
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Table  4 

Test  Configuration  Summary 


Table  5 


Scaled  Total  Fluence  Outputs  of  Various  Laboratory-Scale  TRS  Designs 


Scaled  Total  Fluence 
(cal /gm) 

Mixing 

Technique 

o2/ai 

A1  Particle 

Size  (;j) 

3354 

C 

7.8 

6 

3200 

C 

7.8 

18 

2274 

C 

5.1 

6 

2081 

C 

5.1 

18 

1767 

C 

5.1 

60 

1500 

C 

7.8 

60 
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V SUMMARY  AND  RECOMMENDATIONS 

While  a detailed  theory  concerning  thermochemical  reactions 
(TCRs)  is  not  available  at  the  present  time,  the  basic  laws  of  physics 
seem  to  be  sufficient  to  explain  the  gross  observables  such  as  time  to 
fireball  stabilization,  stabilization  radius,  expansion  energies  and 
energy  radiated  as  a thermal  pulse.  For  many  applications  of  TCRs, 
this  level  of  information  is  sufficient. 

For  the  present  problem,  however,  more  detail  is  necessary. 

The  time  history  of  the  thermal  pulse  and  its  spectrum  must  be  predict- 
able at  least  to  a point  where  the  variables  of  fuel  type,  fuel- 
oxygen  mix,  initial  fuel  charge  and  fuel  feed  rate  can  be  approximately 
determined  before  a test  series  is  begun.  It  is  felt,  at  this  time, 
that  this  level  of  detail  will  be  achievable  with  a more  refined  theory, 
possibly  including  a computer  model  for  fireball  development,  and  data 
which  should  soon  be  available  if  testing  is  continued. 

At  the  conclusion  of  this  investigation  the  following  statements 
summarize  the  state  of  knowledge  concerning  large  scale  optical  flash 
sources : 

1)  Nuclear  waveforms  are  not  directly  simulated  by  any  other 
source.  Their  unique  two  peaked  shape  (time  spectra)  and 
intensity  of  output  are  responsible  for  this. 

2)  The  three  most  promising  sources  of  optical  flash  appear  to 
be  the  SAI  developed  thermochemical  generator  (although  not 
necessarily  using  aluminum  solely),  a high  intensity  flash 
lamp  bank,  and  shock  heated  air  or  argon. 

3)  In  order  to  create  the  double  peaked  time  spectra  two 
sources,  fired  in  rapid  succession,  will  be  required.  Each 
source  will  have  different  properties  to  properly  simulate 
the  first  and  second  maxima. 
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4)  In  order  to  bring  peak  flux,  fluence  and  rise  time  at  a 
distant  detector  into  the  proper  yield  relationship, con- 
siderable "tailoring"  of  the  sources  is  required.  A first 
order  calculation  of  tht  effect  of  source  size  indicated 
that  the  use  of  reflector  panels  to  augment  a small  source 
may  be  as  effective  as  creating  a large  source.  A more 
detailed  cost/benefit  tradeoff  study  will  be  required. 

Freeing  the  source  size  (to  first  order)  in  design  may  allow 
more  readily  attainable  pulse  parameters. 

5)  Concerning  the  SAI  developed  thermochemi cal  optical  pulse 
generation  technique  itself,  we  have  verified  that  the 
prima-cord  is  both  the  probable  ignition  source,  and  the 
cause  of  most  of  the  shock  output.  Reduction  in  prima-cord 
loading  may  prove  advantageous.  The  chemical  reaction  is 
ignited  above  2400°K  (the  melting  point  of  Al^)  and  burns 
initially  at  3600°K  decaying  to  2400°K  at  about  100  msec. 

At  this  point,  the  fireball  is  well  stabilized  and  optical 
radiation  is  relatively  unimportant  to  growth  dynamics.  The 
energy  budget  calculation  at  this  stage  of  development  shows 
that  of  the  initial  energy,  20  percent  goes  to  fireball  expan- 
sion, 70  percent  to  residual  fireball  sensible  heat,  and  10 
percent  to  radiation. 

5)  The  fireball  density  of  both  dust  and  gas  are  calculated 
as  a function  of  time.  At  stabilization  the  fireball  gas 
density  is  about  10  percent  of  normal  atmospheric  density 
and  the  dust  loading  is  10  percent  of  gas  density. 

For  the  tests  conducted  prior  to  1 March  1977,  the  fireball 
is  always  optically  dense,  but  decreasingly  so  as  it  evolves. 
At  stabilization  it  is  three  mean  paths  in  diameter. 
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The  radiative  cooling  rate  calculated  on  the  basis  of  a 
uniformly  heated  sphere  is  consistent  with  measurement. 

Some  additional  cooling  due  to  convective  effects  seems 
to  be  present. 

7)  The  single  theory  proposed  has  been  found  to  agree  quite 
well  with  the  results  of  a much  larger  data  base  than 
that  provided  by  the  fireball  tests  run  by  SAI  so  far; 
that  of  the  fireballs  produced  by  rocket  propellants. 

The  accompanying  figures  (5  and  6 from  Reference  5) 
reproduce  some  data  on  propellant  fireballs  and  include 
points  for  some  of  the  smaller  SAI  thermochemical 
reactions  (fuel  weights  < 5 pounds  in  Figure  5 

and  < 50  pounds  in  Figure  6).  The  points  are  very 

close  to  the  line  extrapolated  from  previous  data  in 

the  first  case  while  in  the  second  the  slope  of  the 

line  for  the  thermochemical  reactions  is  about  the  same 

as  that  for  the  rocket  fuel  case.  Since  the  rocket 

fuel  stabilizes  at  a lower  temperature  than  the  aluminum 

used  in  the  SAI  tests  the  intercept  differences  are  expected. 

8)  At  this  time  it  appears  possible  to  do  some  thermal  pulse 
shaping  by  controlled  injection  of  a metal  powder  fuel 
into  the  fireball  after  initiation.  The  paucity  of  data 
on  such  reactions  has  precluded  a theoretical  prediction 
of  the  effects  such  a procedure  might  produce  however. 

It  would  therefore  appear  to  be  advantageous  to  design 
and  test  such  a device. 

9)  Scaling  laws  evolved  in  the  text  hold  to  within  a factor 
of  two  over  a renion  0.5  to  20,009  gm  of  aluminum.  The 
pulse  duration  for  an  explosively  innited  device  is  pre- 
dictable from  the  scalinn  laws.  Durations  of  the  order 
of  100  msec  seem  feasible.  Additional  scaling  laws  based 
on  ignition  of  premixed  aluminum-oxygen  systems  indicate 
higher  radiant  efficiencies  and  longer  burnino  times  are 
probable.  Considerable  tailoring  of  pulse  duration  seems 
possible. 
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10)  Continued  work  is  required  to  finalize  these  conclusions, 
as  is  noted  above  in  the  area  of  improved  nuclear  wave- 
form and  triggering  requirements  for  NBDS.  More  detailed 
review  of  the  time  history  of  the  multiple  thermochemical 
tests  data  is  required,  and  a thorough  review  of  the 
small  scale  test  is  warranted.  Completion  of  a detailed 
source  augmentation  cost/benefit  study  should  be  accom- 
plished. A detailed  thermochemical  reaction  calculation 
including  collective  effects  should  be  performed.  Other 
particle  types  (magnesium  is  a prime  candidate)  ought  to 
be  considered. 

11)  The  development  of  an  effective  and  inexpensive  NBDS 
simulation  is  now  within  the  range  of  technical  feasi- 
bility. Design  equations  and  scaling  laws  have  been 
developed  in  this  report  for  thermochemical  reactions 
based  on  aluminum  and  oxygen  fuels.  These  coupled  with 
the  NBDS  trigger  requirements  allow  prediction  of  thermo- 
chemical optical  pulse  generators  parameters  that  simu- 
late nuclear  weapon  output.  A detailed  study  of  the 
thermochemical  reaction  process  leading  to  optical  pulse 
output  is  given  in  the  appendix. 
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APPENDIX  A 


1 


THEORETICAL  ANALYSIS  OF  THERMOCHEMICAL  FLASH  GENERATION 


Vredi ction  c f the  behavior  of  -the  thermo chemical  generator 
will  be  based  on  a detailed  understanding  of,  the  physical  and  chemical 
processes  -that  make  up  fireball  development.  Particle  Ignition  and 
burning  -late A may  con.tn.cl  fireball  development  and  radiation  efficiency. 
Fireball  growth  consumer  energy  and  allom  increased  radiation.  The 
details  of  individual  panticle  interaction  help  to  control  output, 
power  and  duration. 


Overall  Models 


In  general,  the  theoretical  knowledge  to  completely  specify 
the  phenomenon  of  thermochemical  flash  generation  via  the  SAI-developed 
optical  flash  unit  is  imperfect.  Considerable  data  exist,  but 
there  is  not  yet  a coherent  theoretical  base  to  support  future  devel- 
opment of  specific  tailored  pulse  output  devices.  It  is  the  purpose 
of  this  section  to  describe  theoretical  models  of  thermochemical  flash 
generators.  Two  models  will  be  developed.  The  first  is  a general 
model  of  fireball  development  that  will  be  useful  for  scaling.  The 
second  type  treats  each  specific  phase  of  the  fireball  development  in 
detail  and  forms  the  basis  of  a detailed  model  of  the  fireball  thermo- 
dynamic development.  The  second  model  may  be  useful  for  detailed 
tailoring  of  the  output  features  of  flash  generators. 


The  fundamental  chemical  reaction  of  burning  aluminum  is: 


2 A1  + | 02  — A1  203  + 7400  cal/gm  al 


Some  of  the  physical  properties  of  Al  and  Al^O^  are  contained 
in  Table  A-l. 


Table  A- 1 


PHYSICAL  PROPERTIES 


AI 

A12°3 


Temperature 
Melting  Boi 1 i ng 

932°K  2740°K 

2318°K  3800 °K 


Latent  Heat 


Fusion  Vapori zation  Density 

94  cal/gm  2619  cal/gm  2.70 

244  cal/gm  --  3.99 


Single  aluminum  particle  ignition  temperatures  in  air  are  measured  to 
be  about  2300°K,  approximately  the  melting  point  of  alumina  (Al^O^). 
Some  controversy  exists  as  to  whether  the  reaction  will  proceed  at 
lower  temperatures  in  dust  clouds  due  to  cooperative  effects. 

The  general  fireball  model  starts  from  basic  physical  princi- 
ples to  develop  the  thermodynamic  relations  for  fireball  size,  dura- 
tion and  growth.  Assume  that  W kilograms  of  fuel  are  available,  and 
at  t = 0 are  ignited.  The  thermochemical  reaction  is  initiated  at  the 
melting  point  of  Al^O^  (^2300°K),  when  the  elemental  aluminum  is  free 
to  react  with  the  oxygen,  and  rapidly  progresses  to  the  vaporization 
temperature  of  Al^O^  (%3800°K).  The  reaction  products  are  assumed  to 
be  gaseous  regardless  of  the  initial  state  of  the  fuel  and  to  remain 
so,  obeying  the  ideal  gas  law,  until  cloud  stabilization.  In  an  uncon- 
fined system,  we  can  take  the  pressure  on  the  gas  to  be  constant  at  pQ, 
atmospheric  pressure,  so  that  the  gas  law  takes  the  form  of 

Vf  . Vi 


NfTf 


N.T. 
1 i 


or 


Nf  T^ 

= HT  TT  - 


N 


f 


(A) 


where  is  the  molar  ratio  of  the  gaseous  product  to  the  gaseous  re- 
actant (i.e.,  1 mole  of  Al^  to  1.5  moles  of  0^)  and  is  the  ratio 
of  the  final  gas  temperature  to  the  initial  gas  temperature  (290°K). 
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The  initial  gas  volume  can  be  determined  by  noting  that  0.88  grams  of 
0^  combine  with  each  gram  of  A1  to  produce  Al^O^,  and  that  at  standard 
temperature  and  pressure  32  grams  of  0^,  in  a pure  gas  form,  has  a 
vol ume  of  22.4  1 i ters . 


We  then  find  that,  in  MKS  units, 
V.  = .616  W 

where  W is  the  weight  of  the  aluminum. 

So  that 


and 


Vf  = 1 .42  x 10  J Tf  W 


rf  = .070  Tf1/3  W1/3  = C Tf1/3  W1/3 


(B) 


= .937  W1/3  at  Tf  = 2400°K 

Support  for  this  simple  model  comes  from  High,  who  reviewed 
data  for  about  60  fireball  radii  produced  by  rocket  propellant  igni- 
tions over  a range  from  1 Kg  to  almost  10®  Kg.  He  found* 

rf(m)  = 1.16  W0-32at  T = 2400°K 

(see  Figure  A-l  taken  from  Reference  7).  The  constant  of  proportional ity 
is,  of  course,  fuel  dependent  which  may  explain  the  difference  between 
the  two  results. 

The  simple  model  can  also  predict  fireball  durations  if  we 
assume  that  the  optical  output  versus  time  is  triangular,  and  define 
the  duration  as  the  time  between  half  power  points.  Then  the  relation- 
ship between  optical  output  (Q),  peak  radiant  flux  per  unit  area  (Q), 
area  (A),  and  time  duration  (t)  is: 

Q = QAT  = dW  (C) 

where  d is  the  optical  output  per  unit  weight  of  A1  (in  cal/Kg). 

*Within  the  error  of  the  figure  and  when  expressed  in  MKS  units 


43 


Peak  flux  per  unit  area  (assuming  radiation  dominates)  is  given  by: 


where 


q = oEd^-T^) 


o = 1.36  x 10  3 cal/m2  sec  °K 


and  E = fireball  emissivity. 

Then  from  Equations  C and  B,  we  find 


" 4 aEcV-T, V/3 

C being  .070  from  Equation  B,  and  that 

t = dW1/3  (F) 

In  evaluation,  t,  C,  and  possibly  E may  be  functions  of  Tp  in  addi- 
tion to  the  explicit  Tp  dependency  shown.  This  in  essence  indicates 
that  t,  the  fireball  duration,  will  be  a function  of  the  cube  root 
of  the  weight  of  fuel.  Reference  5 again  provides  data  on  time 
duration.  Using  the  value  of  2400°K  for  T,  as  a suitable  average 
for  the  test  conditions  under  which  data  was  taken,  and  a value  for 
i of  10^  cal/Kg  and  E = 1 as  representative  of  the  fireball  emissivity 
for  these  fuels  and  C as  indicated  above,  we  thus  obtain 

t(sec)  = 0.201  W(Kg)  1/3  (G) 

while  High  finds  that,  for  rocket  fuels  with  Tp  = 1800  °K 

t(sec)  = 0.180  W(Kg)  0-32 

for  a series  of  about  60  fireball  tests.  Figure  A-2  illustrates  his 
results  from  empirical  data  fitting.  Table  A-2  gives  values  of  the 
coefficient  d and  the  final  temperature  for  the  aluminum  fuel  case. 

The  final  stage  of  development  of  the  model  is  to  convert 
from  a static  to  a time  dependent  or  growth  model.  This  is  accom- 
plished by  dimensional  analysis.  The  expression  for  the  time  dura- 
tion, Eq.  G is  also  true  for  any  instant  of  time,  where  W is  changed 


Table  A-2 

Evaluation  of 
t(sec)  = d W(Kg)  1/3 
As  a Function  of  Final  Temperature 


d 

Tf  (°K) 

0.0302 

3600 

0.201 

2400 

0.470 

2000 

0.769 

1800 
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to  the  fuel  consumed  up  to  time  t.  The  general  equation  for  fireball 
radius  growth  is  then  a function  of  fuel  weight  and  time 

r = f(W,t) 

If  we  assume 

r = fWat'- 

and  use  equations  B and  E and  keep  T constant  then 

CW1/3  = ydM  l+p/3) 
implying  a + 3/3  = 1/3 

While  this  defines  a set  {a,  3}  for  which  the  equation  may  be  valid, 
the  simplest  member  of  the  set  with  non-zero  positive  values  will  be 
taken,  so  that 

a =3=1/4 

and  r = ^(Wt)1^  (H) 

with  Y = C/d 

No  data  from  reactions  other  than  the  Al^Og  system  were  available 
to  test  this  model.  The  agreement  between  the  growth  model  values 

rf(m)  = 2.6  (W(Kg)t(sec) ) 1/4,  Tf  = 3600°  K 

rf(m)  = 1.4  (W(Kg)t(sec))1/4,  Tf  = 2400°  K 
and  the  data  is  acceptable. 

A. 2 Detailed  Fireball  Development  Model 

Basically,  the  Flashbulb  reaction  has  two  distinct  phases: 
the  prima-cord  detonation  phase,  and  the  A 1 /O^  reaction  phase.  At 
t = 0 the  prima-cord  is  ignited,  sending  a shock  wave  through  the 
aluminum  powder  tube,  the  0^  bag  and  into  air  in  the  matter  of  about 
1 msec.  The  shock  heated  and  dispersed  powder  is  accelerated  in  the 
outward  flow  field  and  immediately  encounters  shock  heated  O^.  The 
outer  edge  of  the  expanding  dust  layer  is  ignited  by  this.  The  core 
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area  may  ignite  also  due  to  the  residual  hot  products  of  the  prima- 
cord.  In  any  case  the  dust  layer  is  optically  thick  and  the  surface 
radiation  is  as  from  a black  body.  By  t = 1 msec  the  glowing  dust 
layer  is  visible  in  high  speed  photography.  The  layer  expands  due  to 
drag  in  the  outward  flow  of  shocked  air,  and  due  to  thermodynamic 
expansion.  Radial  stabll ization  occurs  because  of  the  consumption  of  0^ 
within  the  fireball.  Specific  calculations  of  the  following  effects 
were  attempted: 

• Prima-cord  shock  trajectory 

• Prima-cord  shock  velocity 

• Shock  overpressure 

t Shocked  air  temperature 

• Fireball  trajectory 

• Fireball  density 

• Optical  thickness  in  the  Fireball 

• Thermodynamic  expansion 

• Radial  stabilization,  cooling 


A. 3 Shock  Model 


For  the  shock  trajectory  we  used  the  results  of  Jones,  Goyer 
and  Plooster®,  who  show  that  for  a 25  grain  prima-cord  the  shock 
trajectory,  t(r),  can  be  written  as 

t(r)  = 1.27  [(1  + 5.67  x 10-4  r2)1/2  - 1]  msec 
where  r in  cm  is  the  radius  of  the  shock  front  at  time  t. 


Jones  et  al.,  give  a formula  for  the  overpressure  correct  in 
the  weak  shock  limit  as: 

7.7 


Ap 


(1  + 5.8  • 10'4r2)3/8-l 


PS1 


We  also  used  the  approximate  formula 
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for  the  pressure  ratio  across  a strong  shock  in  air,  when  M,  is  the 
flach  number  in  the  undisturbed  region  ahead  of  the  shock.  Generally, 
these  agreed  fairly  well. 

The  temperature  ratio  of  shocked  to  undisturbed  air  was 
estimated  as 


This  is  strictly  valid  only  for  strong  shocks  (Mj  large)  and  indicates 
that  the  shock  wave  is  sufficient  to  ignite  the  A1 ,02  mixture. 

A. 4 Fireball  Parameters  (radius  as  a function  of  time) 

The  fireball  trajectory  was  obtained  empirically  by  measure- 
ment from  high  speed  photography  of  several  Flashbulb  shots.  As 
explained  in  A.l  , these  agreed  with  the  simple  gas  expansion  model. 
Fireball  density  was  calculated  assuming  isotropic  dust  distribu- 
tion. Reference  1 data  for  a variety  of  shots  allows  a formulation 
of  radius  of  fireball  vs.  time  to  be  made  and  is  used  directly  here. 
From  this,  fireball  dust  density  and  optical  thickness  can  be 
found. 

Fireball  gas  density  is  estimated  from  the  perfect  gas  laws 
given  the  measured  fireball  temperature. 

The  comparison  of  the  detailed  quantitative  results  of  Refer- 
ence 1 will  now  be  made.  The  first  general  area  to  be  covered  will  be 
radiative  energy  transfer.  The  relationship  for  fireball  stabilized 
size 

R0(m)  = 2.5  (W(kg) ) 1/3  (J) 

and  fireball  growth 


* 
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R(m)  = 5 (W(kg)t(sec) ) 


(K) 


1/4 

as  determined  in  Reference  1 are  found  to  be  in  agreement  with  the 
overall  model  for  fireball  development  derived  in  Section  (A.l). 

Given  the  relationship  of  fireball  radius  to  aluminum  weight 
and  time  after  ignition,  the  fireball  density  and  optical  thick- 
ness will  now  be  determined.  The  above  equations  imply  that  the  time 
of  fireball  radial  stabilization  is 

tQ(seconds)  = 0.0625(W(kg) )1/3  (L) 

and  the  density  of  the  fireball  gas,  fireball  dust  and  the  fireball 
optical  thickness  can  be  calculated. 

The  fireball  gas  density  is  given  by  the  ideal  gas  law  - using 
the  stabilization  temperature  of  2400°K 

Pgas  = .157  mg/cc,  at  fireball  radial  stability 
Pgas  = .100  mg/cc  at  3600°K  at  the  peak  power  point 

The  average  dust  density  at  fireball  radial  stabilization  (assuming  homo- 
geneous mixing)  is  given  by 

Mass  of  Aluminum 
' dust  Volume  of  Fireball 

pdust  = = "~4 — KS'  3 = 0.015  mg/cc  at  stabilization 

3^0  ^(2.5)  nr 

p , ,/p  = 0.1  at  stabilization 

dust  gas 

The  dust  density  as  a function  of  time  can  be  found  by  application  of  the 
fireball  growth  equation  rather  than  using  the  radius  at  stabilization. 

The  more  general  equation  is  then 
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_ 0.0019  ju(kg)l  1/4  /mt[\ 

t(  } ' (t(sec))3/4  W 


' dusl 


The  optical  thickness  is  calculated  from  the  number  density  per  cubic  cm., 
N,  and  the  optical  cross  section  o,  thus 


l 


optical  No 


using  the  number  density  definition, the  average  dust  density  is  defined 
as, 

pdust=  Nxpalx  (^3) 

when  a is  the  average  particle  radius,  and  pfll  is  the  density  of  solid 
aluminum.  Thus, 

'optical  Pdusto 

where  a is  the  particle  radius  in  cm. 

2.7  gm  Y 4 33  A 

_3  X cm  x 104  p/cm 

it. 


cm 


'optical  31/4  „ „ 

^0019x10  V'A  (M  \ 2,a2c»l 2 

t3'4  ' cm3/ 


£(m)  =0.94  a(p ) 


t(sec) 


3/4 


W(kg) 


-1/4 


l =0.94  m for  a = Ip,  t = lsec. 


= 5.70  m for  a = 6p,  t - lsec 


The  following  table  (Table  4-3)  indicates  the  values  of  fireball 
;i ze,  dust  density  and  optical  thickness  as  a function  of  time  for  1 Kg, 
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aluminum  dust  of  3 micron  radius  particle  size.  The  cloud  is  assumed  to 
be  at  3600°K  to  after  10  3 seconds,  and  2300°K  for  10~3  seconds  and 
later.  Reference  to  the  table  will  show  that  the  fireball  is  always 
optically  thick.  For  example,  it  is  12  mean  f t ee  paths  thick  at  0.1 
seconds  well  after  all  significant  optical  output  has  ceased.  An 
expression  for  the  ratio  of  the  fireball  radius  to  optical  thickness 
can  be  derived  from  the  general  formula  above.  The  optical  thickness 
at  stability  is 

£0(ni)  = 0.12  a (p ) 
independent  of  weight  of  aluminum. 

The  ratio  of  fireball  radius  to  optical  thickness  at  stabili- 
zation is  given  by 

= 21 . 3(W(kg ) ) 1/3 

a0  3 (p  ) 

Table  A-4  shows  the  results  of  the  evaluation  of  this  formula  for  a variety 
of  particle  sizes  and  aluminum  weights.  These  tables  indicate  that 
under  the  test  conditions  of  Reference  1 ( a = 3y  W >_  500gm  ) the 
fireball  is  of  sufficient  optical  thickness  for  good  radiation  efficiency. 

A. 5 Fireball  Expansion 

The  second  major  energy  loss  mechanism  to  be  studied  is  that  of 
fireball  expansion.  The  fireball  expands  because  the  rate  of  energy 
generation  exceeds  the  ability  of  the  surface  to  radiate.  Thus  the  fire- 
ball expands,  absorbing  energy  in  thermodynamic  work,  and  increasing  the 
surface  to  allow  more  radiation.  Stabilization  is  reached  when  the  rate 
of  generation  and  loss  equalize.  The  fireball  radius  is  given  by  equa- 
tions J and  K.  The  energy  absorbed  in  expansion  after  initiation  is 
given,  in  general , by 

Ee(t)  = pQ  (Vf(t)  - (M) 

1 P0  vf(t) 
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where  pQ  = local  atmospheric  pressure  and  is  the  fireball  volume  at 
time  t.  Thus  the  total  energy  consumed  in  expansion  to  radial  stabiliza- 
tion volume  VQ  is 

Eeo  = poVo  = L57  x 1q3  cal/9mA£  00 

The  expansion  energy  as  a function  of  time  is  given  by 

Ee(cal)  = 1.26  x 107  W(kg)t(sec)  3/4  (0) 

and  the  power  absorbed  in  expansion  is  given  by 


Pe(cal/sec)  = 


9.45  x Hr  W( 
t(sec)  1/4 


while  the  power  loss  per  unit  area  is  given  as 


P (cal/cni  sec) 


3.01  (U(kg))1/4 


(t(sec) 


The  power  loss  due  to  expansion  is  given  in  Table  A-5  for  various 
times  and  aluminum  weights.  Since  pr  the  radiative  power  per 
unit  area  (assuming  a black  body  and  unit  emissivity)  is  given  by 

n -r-4 

Pr  = oT  (R) 

with  a = 1.36  x 10  ^ cal/cm<Lsec°K4 

the  growth  time  to  the  point  when  radiation  and  expansion  losses  (R  and  Q) 
equalize  is 


t =/Ml 

B U4 


( W ( kg ) ) 
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A final  calculation  can  be  made  of  the  rate  of  consumption  of 
the  aluminum  fuel.  The  output  rate  is  made  up  of  radiative  and  expan- 
sion losses  together  must  equal  that  of  fuel  consumption.  From 


where  u is  the  effective  value  of  the  fuel  (^3 106  cal/kg)  and  tr  is  the 
net  consumption  time  constant.  Then,  using  Eqs.  P and  D with  the 
area  of  the  fireball  surface,  we  have 


4TOT4(5W1/4t1/4)2 


+ 9.45  x 106  - 
t 


3/4 

174 


~ = 231  (^) 1/2  + 3.04(^)1/4  at  T = 3600°K 


Table  A-6  lists  values  of  t at  peak  output  for  various  aluminum  weights. 
Also  shown  is  the  particle  radius  a that  has  the  same  ignition  and 
burning  time,  and  the  rate  of  consumption  of  the  fuel  W = W/t  . 


Table  A-6 

Results  of  Aluminum  Consumption  Calculation 


w 

tr 

a 

W 

kg 

sec 

p 

kg/sec 

0.001 

0.004 

13 

0.25 

0.01 

0.008 

22 

1.25 

0.1 

0.018 

35 

5.60 

1.0 

0.C40 

55 

25.00 

10.0 

0.080 

80 

125.00 

This  particle  radius 

is  therefore 

di fferent 

from  that  shown  in  Table  3 

which  was  determined 

by  the  time 

to  fireball 

peak  output  rather  than  the 

longer  time  for  all  radiation  output  used  here. 
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The  thermodynamic  expansion  of  the  fireball  is  calculated 
via  the  formula 


“th  ■/ 


pdv  = P . (V,  - V.) 

r atmos  f 1 ' 


V. 

1 


for  the  case  of  expansion  of  a gas  cloud  in  the  atmosphere. 

The  cooling  rate  of  the  fireball  is  determined  by  the 
equation  for  the  rate  of  energy  loss: 


4 3 dT  _ i-t-4  . 2 

7 n r pc  -nr=cET  4Tir 
3 p dt 


as  shown  in  equation  (3),  page  14,  of  Reference  2 which  yields 

dT  = 3c  ET4 
dt 


(T) 


rpcp 

where  p and  c^  are  the  density  and  the  specific  heat  at  constant 
pressure  for  the  fireball.  The  cooling  rate  of  a spherical  shell 
would  be  of  similar  form  but  with  no  1/r  dependency. 

The  final  calculation  that  will  be  made  here  is  of  the  cooling 
rate  of  the  fireball  using  Eq  T. 

r4 


dT  _ 3g  E T 
dt 


usi ng : 


rpc 


cr  = 1.36  x 10' 


•12 


cal 


2 o,. 

cm  sec  K 


E = 1 

Cp  = 6 cal/mole-°K 


and 

we  find 


pfi reball  = 0.4g/ltr  at  2300°K 
^ = 15000  ~ at  T = 2300°K 


A value  of  about  20000°K/sec  at  2300°K  for  fireball  cooling  was  observed. 
Convective  or  turbulent  heating  of  the  air,  test  objects  and  the  ground 
surface  could  account  for  the  additional  cooling. 
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A. 6 Energy  Balance 

Calculation  via  energy  balance  of  the  total  energy  radiated 
per  gram  of  aluminum  consumed,  up  until  the  stabilization  radius  is 
achieved  at  2400°K,  can  now  be  made.  From  conservation  of  energy 
we  find  that  E-j  the  energy  liberated  is 

E]  = ^ai  ^a I -'T  + ^air  ^air  AT  + La]  Mai  - 

+ E + E 
e r 

where  the  first  two  terms  result  from  the  specific  heat  contributions 
of  the  aluminum  and  free  air  in  the  fireball,  the  next  two  result 
from  the  latent  heat  of  fusion  of  the  aluminum  and  A^O^  (the  last 
latent  heat  is  liberated  as  the  A1  ^0^  fuses  and  hence  is  negative) 
and  the  last  two  terms  represent  the  energy  lost  in  expansion  of  the 
fireball  and  in  radiation. 

The  first  term  is  simply 

C , M.  AT  = .27  cal/gm°Kxlgm  x 2110°K  = 570  cai 

a I a I 

For  the  second,  we  note  that  the  mass  of  air  at  2400°K 
fills  a volume  of 

V2  = 4/d  n R03 

while  from  the  ideal  gas  law  its  original  volume  was 

Q 

V1  = V2  Tl/T2  = 4/3  11  Ro3  Iwo0  = 8‘54  x 10"3m3 

or  .353  moles  of  air  so  that,  using  a molar  specific  heat  we  find 

C,,-,,  AT  = 6 cal/mole°K  x .353  moles  x 2110°k  = 4470  cal. 
dir  d»r 

The  third  term  is  just  the  latent  heat  of  fusion  of  one  gram  of  alum- 
inum or  94cal  and  for  the  fourth  we  use  the  latent  heat  of  fusion  of 
A1  ,0^  (244  cal/gm)  and  the  fact  that  one  gram  of  A1  produces 

1.88  gm  Alo0  so  that  this  term  is 
3 2 3 

- Lai  „ M, i n = -244  cal/gm  x 1.88  gm  = -4b9  cal 
a 2°3  “23 
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The  fifth  term  is  obtained  from  the  previous  equation  for  Eg 
Ee  = 1.26  x 107  [Wt]3/4 

with  W = 10  3 Kg  and  t = t the  stabilization  time  = .0625  W373  so  that 

E = 1575  cal 
e 

Thus  using  the  previous  value  of  E^  (7400  cal/gm)  we  find 
7400  = 570  + 4470+  94  - 459  + 1575  + ER 
or  Er  = 1150  cal/gm 

Considering  that  this  approach  has  used  average  values  for  the 
temperature  (and  phase)  dependent  specific  heats  and  has  assumed  any 
processes  occurring  above  24Q0°K  are  reversed  without  loss  of  energy, 
this  value  is  in  very  good  agreement  with  the  observed  value  of  approx- 
imagely  1100  cal/gm. 

A. 7 Dust  Explosion  or  Fire 

The  final  areas  of  interest  are  those  of  the  burning  rate 
temperature  and  dust  flame  propagation  in  aluminum/oxygen  mixtures. 
References  9and  10  describe  a series  of  experiments  on  aluminum  particle 
thermodynamics.  This  discussion  will  cover  three  aspects;  single 
particle  ignition  and  burning,  particle  accretion,  and  dust  ignition 
both  in  deflagration  and  detonation  modes. 

The  ignition  temperature  of  the  single  aluminum  particle  depends 
on  the  partial  pressure  of  oxygen  varying  from  about  2300°K  at  air  at 
1 atmosphere  (P02  = 0.16)  to  about  2150°K  at  P02  = 1.0.  Ignition  temper- 
atures of  single  carbon,  and  magnesium  particles,  vary  as  1/r  but  aluminum 

does  not  appear  to  be  r dependent.  Due  to  general  heat  transfer  arguments, 

9 

ignition  times  for  all  particles  should  vary  as  r . Burning  life-times 
vary  as  r due  to  the  cenosphere  type  of  burning  (diffusion  through 
A I2O3  shell)  and  vary  inversely  with  the  oxygen  diffusion  coefficient. 

This  in  turn  is  linearly  proportional  to  oxygen  partial  pressure. 
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Figure  A-3  (taken  from  reference  9)  shows  data  on  aluminum 
particle  ignition  times.  From  this  data  a general  equation  for  the 
ignition  time  can  be  extracted: 

t.(sec)  = 2.5  x io'6  (a(M))2  (U) 

Burning  times  in  infinite  excess  oxygen  can  be  estimated  from  limited 
data  as 

tg(sec)  ^[0.0002  a (p ) - .0004]  (V) 

Burntime  approximately  doubles  at  50.  excess  oxygen.  Table  A-7  gives 
some  data  from  the  literature. 

Although  not  completely  verified,  the  general  features  of 
cenospheric  particle  burning  appear  to  be  valid  for  the  case  of  aluminum 
in  air.  Figure  A-4  (taken  from  reference  9)  illustrates  the  features  of 
interest.  Diffusion  of  gaseous  aluminum  to  the  combustion  zone  through 
the  porous  alumina  shell  proceeds  until  the  aluminum  is  consumed. 

Figure  A-5  (taken  from  reference  9)  plots  the  ignition  and  burning  condi- 
tions as  a function  of  partial  pressure. 

When  a collection  of  particles  (a  dust  cloud)  is  considered  there 
are  cooperative  effects  that  may  take  place.  The  ignition  temperature 
of  dust  clouds  generally  varies  as  1/p  at  least  in  the  range  of  1/10 
air  density.  This  may  be  due  to  several  effects:  Brownian  motion, 
turbulent  density  fluctuations,  and  accretion.  The  overall  tendency 
would  be  to  produce  larger  particles  and  of  higher  local  density  which 
may  be  easier  to  ignite. 

Particle  burning  lifetime  is  extended  in  a dust  cloud.  While 
theory  shows  that  an  exact  stoichiometric  ratio  yields  an  infinite 
burning  time,  a stoichiometric  ratio  of  1.1  (oxygen  excess)  lengthens 
burning  by  nearly  4 times  over  a single  particle  case.  During  burning 
the  oxygen  partial  pressure  is  reduced  and  this  effect  lengthens  the 


Figure  A-3.  Experimen 
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The  formulae  describing  all  these  phenomenon  are  complicated  but 
tractable.  They  will  not  be  repeated  here  but  are  given  in  References 
8 , 9 and  10. 

Dust  cloud  ignition  leads  to  deflagration  (flame)  or  detonation 
(shock).  In  deflagration  the  burning  velocities  are  of  the  order  of  1 
up  to  1,000  m/sec,  while  detonation  proceeds  at  2000  to  6000  m/sec,  and 
are  accompanied  by  shocks.  The  various  dependencies  described  are 
tabulated  in  TableA-8.  Various  types  of  particles  may  be  considered. 

For  example,  magnesium  burns  in  about  one-half  the  time  of  aluminum. 

Both  carbon  and  magnesium  have  considerably  lower  ignition  temperatures . 
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Table  A-8.  Parametric  Variations  of  Particle  Thermodynamic  Variables 


r 


Single  Particle 

Dust 

Cloud 

Og  Partial  Particle 

Cloud 

Stoic. 

Pressure  Radius 

Density 

Ratio 

P 

0 

r 

6 

s 

Ignition 
Temp.  T. 

L 

P"2 

0 

r'1 

A 

P-1 

N/A 

Igni tion 

Time  t^ 

w 

2 

r 

N/A 

N/A 

Burn 

Time  t^ 

Po! 

2 

r 

N/A 

-Ins 

A)  Verified  for  carbon  and  magnesium. 

B)  Depends  on  thermal  conductivity  of  gas  and  . 
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